On the Cooling of the Neutron Star in Cassiopeia A 
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We demonstrate that the high-quality cooling data observed for the young neutron star in the 
supernova remnant Cassiopeia A over the past 10 years-as well as all other reliably known tem- 
perature data of neutron stars-can be comfortably explained within the "nuclear medium cooling" 
scenario. The cooling rates of this scenario account for medium-modified one-pion exchange in dense 
matter and polarization effects in the pair-breaking formations of superfluid neutrons and protons. 
Crucial for the successful description of the observed data is a substantial reduction of the thermal 
conductivity, resulting from a suppression of both the electron and nucleon contributions to it by 
medium effects. We also find that possibly in as little as about ten years of continued observation, 
the data may tell whether or not fast cooling processes are active in this neutron star. 

PACS numbers: 97.60.Jd, 95.30.Cq„ 26.60.-c 



Introduction.- The isolated neutron star in Cas- 
siopeia A(Cas A) was discovered in 1999 by the Chandra 
satellite [IJ. Its association with the historical supernova 
SN 1680 [2| gives Cas A an age of 330 years, in agree- 
ment with the nebula's kinematic age [3j. The distance 
to the stellar remnant is estimated to be 3.4^o 3 kpc Q. 
The thermal soft X-ray spectrum of Cas A can be fitted 
with a non-magnetized carbon atmosphere model, a sur- 
face temperature of 2 x 10 6 K, and an emitting radius of 
8 to 17 km (Ht Analyzing the data from 2000 to 2009, 
Heinke & Ho [6J reported a rapid decrease of Cas A's sur- 
face temperature over this 10-year period, from 2.12 x 10 6 
to 2.04 x 10 6 K. Such a rapid drop in temperature con- 
flicts with standard cooling scenarios based on the effi- 
cient modified Urea (MU) process 0, Q . First interpre- 
tations of Cas A's temperature data were provided very 
recently by Page et al. @ and Yakovlev et al. lfl 11|. 

The interpretation of Page et al. Q is based on the 
"minimal cooling" paradigm [l2j , where a minimal num- 
ber of cooling processes is taken into account. These 
are photon emission, the MU process, nucleon-nucleon 
(NN) bremsstrahhmg (NB) and the neutron (n) and 
proton (p) pair breaking-formation processes (nPBF and 
pPBF). The latter are particularly important in the 
ultra-dense cores of neutron stars ljl-15|, where neu- 
trons form Cooper pairs in the 3 P 2 channel and proton 
pairing occurs in the 1 5o channel. To calculate the NN 
interaction entering the emissivities of the MU and NB 
processes the minimal cooling scenario employs the free 
one-pion exchange (FOPE) model As shown in Q, 
the Cas A data can be neatly reproduced by assuming a 
large value for the proton pairing gap throughout the en- 
tire stellar core and by fixing the critical temperature for 
the neutron 3 P 2 pairing gap at around 0.5 x 10 9 K. The 
result is mildly sensitive to the neutron star mass. Sur- 



face temperature-age data of other neutron stars, which 
do not lie on the cooling curve of Cas A, are explained 
within the minimal cooling scenario mainly by assuming 
variations in the light element mass of the envelopes of 
these stars. 

The work of Yakovlev et al. 1(J 11 | includes all emis- 
sion processes which are part of the minimal cooling 
paradigm and uses also the FOPE to model the NN 
interaction. As in Q, it is assumed that the proton 
gap is large and non-vanishing in the entire stellar core. 
The latter assumption facilitates a strong suppression of 
the emissivity of the MU process. The value and the 
density dependence of the 3 P 2 neutron gap are fitted 
to the Cas A data, leading to a critical temperature of 
(0.7 — 0.9) x 10 9 K for the neutron pairing gap. Both 
groups therefore came to the striking conclusion that the 
temperature data of Cas A allow one to extract the value 
of the 3 P 2 neutron pairing gap. 

In this Letter, we present the "nuclear medium cooling 
scenario" as an alternative model for the successful de- 
scription of the temperature data of Cas A. Aside from 
describing the Cas A data extremely well, this model 
reproduces also all other presently known temperature 
data of NSs, without the need of making any additional 
assumptions. Before representing the stellar cooling re- 
sults, we outline the key features of the nuclear medium 
cooling scenario next. 

Nuclear medium cooling.- Motivated by the fact that 
the existing temperature-age data of neutron stars seem 
to be incompatible with a unique cooling evolution, the 
nuclear medium cooling scenario has been worked out in 
12 "HI- It provides a microscopic justification 



Refs. 14 



for a strong dependence of the main cooling mechanisms 
on the density (and thus on the neutron star mass) . The 
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nuclear medium cooling scenario has been successfully 
applied to the description of the body of known surface 
temperature-age data of neutron stars [H, 0, Hl| . The 
scenario addresses the often disregarded role of medium 
effects on the MU and NB processes. Furthermore, as it 
is commonly accepted, the neutron and proton superflu- 
idity with density dependent pairing gaps is causing an 
exponential suppression of neutrino emissivities of the 
nucleon processes and of the nucleon specific heat, and 
opens up the new class of nPBF and pPBF processes. 
We also want to stress that the thermal conductivity is 
essential for the cooling of young objects such as Cas A. 
The next paragraph is devoted to a brief discussion of 



the medium modified Urea (MMU) to the MU process, 



these issues. For more details, we refer to 181-2 



1. Free versus medium-modified one-pion-exchange in 
dense matter: The insufficiency of the FOPE model for 
the description of the TVTV-interaction is a known is- 
sue [3]. Indeed, calculating the MU emissivity per- 
turbatively one may use both the Born TV TV interaction 
amplitude given by the FOPE and the imaginary part 
of the pion self-energy. In the latter case one needs 
to expand the exact pion Green's function D^{ui,k) = 
[u 2 — m 2 — k 2 — H(u>, fc,™)]" 1 to second order using for 
the polarization function II(u;, k, n) the perturbative one- 
loop diagram, ITo(u;,fc, n). For k = ko, which is the 
pion momentum at the minimum of the effective pion 
gap uj* 2 = — Z? 7i r 1 (cj = 0, k = fc ), the polarization func- 
tion IIo(a;, A - = fco — PF,n, n) yields however a strong TV TV 
attraction. (Here is the pion mass and pf,u is the 
neutron Fermi momentum.) This attraction is so strong 
that it would trigger a pion condensation instability al- 
ready at low baryon densities of n ~ 0.3 no, which is 
in disagreement with experimental data on atomic nu- 
clei, (no = 0.16 fm -3 denotes the nuclear saturation 
density.) The discrepancy is resolved by observing that 
together with a pion softening (i.e., a decrease of the ef- 
fective pion gap u*(n) with increasing density) one needs 
to include the repulsion from the dressed irNN vertices, 



ith C ~ 1.4-=- 1.6. A consis- 



r ~ [1 + C(n/n o yi*]- 
tent description of the TV TV interaction in matter should 
thus use a medium modified onc-pion exchange (MOPE) 
interaction characterised by the full Green function of 
the dressed pion, dressed vertices T(n), and a residual 
TVTV interaction, as done in this Letter. According to 
17 . HH, the main contribution for n > uq is given by 
MOPE whereas the relative contribution of the residual 
interaction decreases with increasing density. Following 
the model used in [20I 21 ] pion condensation may arise 
only for n > n^ r = 3 no, i.e., for NS masses M > 1.32M 
within a rclativistic version of the APR equation of state 
[22I ] which we use. In the calculation of the neutrino 
emissivity not only radiation from the nucleon legs but 
also from intermediate reaction states is now allowed. 
With such an interaction the ratio of the emissivity of 
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strongly increases with density (for n > no). Although an 
increase of the ratio of emissivities of the medium modi- 
fied nucleon (neutron) brcmsstrahlung process (MnB) to 
the unmodified brcmsstrahlung (nB) is less pronounced, 
the MnB process, being not affected by the proton super- 
conductivity, may yield a relatively large contribution in 
the region of a strong proton pairing. 

2. Pair-breaking formation: The important role of po- 
larization effects in pPBF and nPBF processes was first 
noted in [lij ]. Recently, additional support came from 
an analysis of the vector current conservation in PBF 
reactions [3, 0). In these reactions, diagrams with 
the normal and anomalous Green functions turn out to 
cancel each other, so that the main contribution to the 
PBF emissivity comes from processes of the axial current 

24l | . Another important in- medium effect was recently 
observed in the calculation of the neutron pairing gap 
3 P 2 . Taking into account the polarization effects, it was 
shown in [25[ that the associated gap A rm ( 3 P 2 ) ~ keV, 
i.e. it is dramatically suppressed compared to BCS based 
calculations 26[. For completeness, we also mention the 
possibility of a strong enhancement (more than 1 McV) 
of the gap, as argued in 27J. At first glance, the results 
of [2H and (22| seem to illustrate uncertainties in value 
of the 3 P2 gap, which would suggest to treat A ran ( 3 P 2 ) 
as a free parameter in cooling studies. This, however, 
is not the case since the solution of the gap equation of 

27| exists only for A„„( 3 P 2 ) ~ 1 MeV and disappears 
for smaller values of the gap. Moreover, they use the 
approximation < uj* 2 (n) <C m 2 , so that their new so- 
lution may exist only within a narrow range of the criti- 
cal point for the onset of pion condensation. In realistic 
treatments, pion condensation appears always as a first 
order phase transition 18j with a jump of lo* 2 from a 
positive to a negative value. The required small values 
of the pion gap may therefore not be achieved. More- 
over, [21| has verified that the cooling data are hardly 
described if the gap A„„( 3 P 2 ) were large (> 1 MeV) 
over a broad density region. We therefore disregard the 
possibility of a large value of A„ n ( 3 P 2 ) and adopt tiny 
A lm ( 3 P 2 ) following [25j|. The A„„( 1 S , ) proton gap is 
taken from [2^]. Two different models, labeled I and II, 
are used for the A pp ( 1 S'o) gap [13, HH- Model I is from 
[29I ] and model II is from the calculations in [26| . NS cool- 
ing data can be well described within the nuclear medium 
cooling scenario for both models I and II [2(j, provided 
the A nn ( 3 P 2 ) is str ong ly suppressed in agreement with 
Schwenk & Friman [25j] . 

3. The heat conductivity: The heat conductivity, K, of 
supcrfluid neutron star matter is another key ingredient 
crucial for the cooling of young neutron stars, such as 
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Cas A. It is given by k =JCi K i where ft, are the par- 
tial contributions to k. In [5(| and the electron heat 
conductivity n e computed according to [30[ was used. 
More recent studies [3l| showed that for temperature and 
density regions where nucleons are non-superfluid, Lan- 
dau damping reduces K e below n n . Proton superfluidity 
causes a Meissner screening of transverse photons. Tak- 
ing all these effects into account reduces the total ther- 
mal conductivity [3l[ by an order of magnitude. Also, 
in [20I ] we argued that medium effects may significantly 
suppress the neutron contribution K n to the thermal con- 
ductivity, the effect being not included in j3l[. Indeed, 
K n oc 1/|M 2 |, where M is the NN interaction matrix el- 
ement discussed already in the MOPE context as being 
enhanced for n > no. The impact of a low thermal con- 
ductivity on the thermal evolution of neutron stars ac- 
complished by introducing a factor £ K = 0.3 was demon- 
strated in Fig. 17 of [20|. The net effect is a delay of the 
temperature decline of young (~ 300 yr) neutron stars. 
This idea of a possible strong suppression of the ther- 
mal conductivity, as supported by [3l|, proves essential 
for the explanation of the rapid cooling of Cas A in this 
letter. 
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FIG. 1: (Color online) Cooling of NSs with nuclear medium 
effects, with and without pion condensation (PU), see also 
Fig. 17 of [13]. Data from Refs. QUEl. 



The Neutron Star in Cas A.- The ingredients of the 
nuclear medium cooling scenario discussed above lead to 
neutron star cooling curves in Fig. 17 of Ref. [2(| where 
model I for the proton gap has been used and the role of 
the heat conductivity on the hot early stages of hadronic 
neutron star cooling was elucidated. In Fig. Q] we redraw 
those cooling curves allowing for a minor readjustment of 
the heat conductivity parameter. The bold curves arc for 
a heat conductivity suppressed by a factor of £ K = 0.265, 
while the thin lines are for the unsuppressed heat con- 
ductivity of One sees that for a suppression factor 
of ( K = 0.265 and a stellar mass of M — 1.463 M Q (blue 
bold solid line) we arc able to fit the temperature data 
for Cas A perfectly, as can be seen from the magnified 
10-year epoch for which high-precision cooling data ex- 
ist. This star is our best-fit model. Lowering the neutron 
star mass to M = 1.390 M Q (red dash-dotted line) the 
whole set of available cooling data is covered. Assuming 
the absence of a pion condensate in the core of a neutron 
star, the Cas A cooling data can still be reproduced by 
reducing £ K from 0.265 to 0.220 and readjusting the neu- 
tron star mass to a somewhat higher value of 1.506 Mq, 
sec Fig. 1. The proton gap of model II is significantly 
smaller than that of model I. Nevertheless, the Cas A 
data can still be nicely fitted for £ K < 0.015 and neutron 
star masses M > 1.73 Mq. 

To demonstrate the impact of the heat conductivity on 
the cooling process we present in Fig. [2] the temperature 
profiles for the 1.463 Mq neutron star (f K = 0.265) for 
stellar ages from 10~ 8 to 10 3 years. One sees that the 
heat conductivity is important during the first t < 300 
years and would thus affect the cooling history of Cas A. 
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FIG. 2: (Color online) Snapshots of temperature profiles for 
the Cas A cooling curve (blue bold solid line) of Fig. [TJ 

In Fig. [3] we show the individual contributions of the 
cooling processes of our scenario to the total neutron 
star luminosity for the neutron star, M = 1.463 Mq 
and = 0.265, which best reproduces the cooling of 
Cas A in Fig. [TJ We see that the nMMU is the most effi- 
cient process in our scenario, while all PBF processes are 
less important. The MnB and MpB luminosities domi- 
nate over those of PBF. They are not shown in Fig. [3] 
since they have rather similar shapes as the nMMU and 
pMMU curves. Note that PU processes affect the NS 
cooling primarily at later times. 

Summary and Conclusion. -We have shown in this Let- 
ter that the nuclear medium cooling scenario allows one 
to nicely explain the observed rapid cooling of the neu- 
tron star in Cas A. As demonstrated already in [2(J, in 
this scenario the rapid cooling of very young objects like 
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FIG. 3: (Color online) Individual contributions of the cool- 
ing processes nMMU and pMMU, 1S pPBF and nPBF, 3P 2 
nPBF, PU, and surface photon emission to the total stellar 
luminosity for the neutron star shown in Fig. [5] 



Cas A is due to the efficient MMU and MnB processes, 
a very low (almost zero) value of the 3 P2 neutron gap, 
a large proton gap and a small thermal conductivity of 
neutron star matter. 

Our explanation of the Cas A cooling constitutes an 
alternative to that of 0, HH , which is based on a strong 
PBF process due to 3 P2 superfluidity in neutron star in- 
teriors. We support, however, the conclusion of these au- 
thors that a large value of the proton gap is preferable, 
albeit not necessarily in the entire neutron star core. The 
results presented in Fig. Q] predict that the rapid cooling 
observed for Cas A will continue for a few more decades 
until it slows down when the temperature domain around 
log 10 T s [K] = 6 is reached. Already in about ten years 
from now of continued monitoring, the high accuracy of 
the data for Cas A's surface temperature will allow one 
to distinguish at the 2 a level between models with and 
without additional fast cooling processes (pion conden- 
sation in our case). 

To discriminate between alternative cooling scenarios 
further tests may be considered, such as the comparison 
of log N-log S distributions from population synthesis 
with the observed one for isolated neutron stars. A re- 
cent study of this kind [32| favored model II for the pro- 
ton gaps. Thus it may well be that actual values of the 
thermal conductivity are smaller than assumed in Fig. 1, 
or that there are other important aspects of the cooling 
of Cas A which have not yet been identified. 
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